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ABSTRACT 

The  fast  and  slow  cookoff  behaviour  of  two  series  of  candidate  insensitive 
booster  compositions  based  on  RDX/Elvax  210,  and  incorporating  various 
amounts  of  PETN  and  TATB,  has  been  numerically  simulated  using  a  one¬ 
dimensional  finite  difference  code.  The  model  solves  a  cylindrically  symmetric 
heat  flow  equation  for  a  mixture  of  two  energetic  materials  with  a  time 
dependent  boundary  temperature.  The  temperature  dependence  of  the  thermal 
conductivity  and  specific  heat  of  each  of  the  explosives  is  included,  as  well  as  the 
effect  of  melting,  and  the  effect  of  different  kinetic  schemes  for  the 
decomposition  of  the  RDX.  The  simulations  accurately  reproduce  the  time  to 
explosion  and  surface  temperature  at  explosion  for  varying  PETN  concentration 
at  both  fast  and  slow  heating  rates,  and  also  provide  good  agreement  with 
experiment  for  varying  TATB  levels  at  the  slow  heating  rate.  However  at  the  fast 
heating  rate  there  is  a  divergence  between  the  simulated  results  and  experiment. 
The  calculations  clearly  illustrate  the  need  to  include  the  temperature 
dependence  of  the  thermal  properties  of  the  material,  and  a  kinetic 
decomposition  scheme  appropriate  to  the  degree  of  confinement,  before  good 
agreement  between  simulated  and  experimental  results  can  be  obtained. 
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Simulation  of  Cookoff  Results  in  a 
Small  Scale  Test 


EXECUTIVE  SUMMARY 


The  adoption  of  an  Insensitive  Munitions  policy  by  the  ADF  has  led  to  the 
requirement  to  develop  an  insensitive  booster  composition.  A  series  of 
compositions  based  on  the  common  military  explosive  RDX,  coated  with  an 
ethylene-vinyl  acetate  copolymer  Elvax  210,  and  containing  varying  amounts  of 
the  explosives  PETN  or  TATB,  has  been  tested  experimentally  for  shock,  impact 
and  thermal  sensitivity.  None  had  sufficiently  low  thermal  sensitivity  to  qualify 
as  an  insensitive  booster  composition. 

To  understand  the  thermal  behaviour  of  these  materials  in  more  detail  a 
mathematical  model  of  the  heat  flow  within  the  candidate  booster  compositions 
was  developed  and  implemented  in  a  one-dimensional  computer  code.  The  code 
was  then  used  to  simulate  the  experimental  test  procedures.  It  was  found  that 
the  code  could  accurately  predict  the  heat  flow  within  the  compositions  in  most 
cases,  provided  that  allowance  was  made  for  the  temperature  dependence  of  the 
thermal  properties  of  the  explosives.  It  was  also  found  that  the  rate  of 
decomposition  of  the  explosives  depended  strongly  on  the  degree  of 
confinement,  and  that  a  kinetic  decomposition  scheme  appropriate  to  the 
strength  of  the  confinement  was  needed  before  good  agreement  with  experiment 
was  obtained.  The  material  models  developed  in  this  study  can  now  be  used  in 
more  advanced  computer  codes  which  attempt  to  predict  the  degree  of  violence 
of  the  thermal  response. 
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1.  Introduction 


A  recent  MRL  Report  described  the  preparation  and  characterisation  of  three  series  of 
candidate  insensitive  booster  compositions  based  on  RDX/Elvax  210  and 
incorporating  various  amounts  of  PETN,  DATB  or  TATB  [1].  Each  composition  was 
characterised  for  impact  sensitiveness,  shock  sensitivity  and  cookoff  behaviour.  One 
of  the  compositions  had  suitable  shock  and  impact  sensitiveness,  but  none  of  them 
exhibited  sufficient  reduction  in  cookoff  reaction  violence  to  give  an  acceptable 
insensitive  booster  composition. 

In  an  effort  to  gain  a  greater  understanding  of  the  behaviour  of  energetic  materials 
when  subjected  to  heating  at  various  rates,  we  have  numerically  simulated  the  thermal 
environment  in  the  Super  Small-scale  Bomb  Test  (SSCB)  for  two  of  the  candidate 
insensitive  booster  compositions  described  by  Dagley  et  al.  [1].  Our  model  solves  the 
one- dimensional  cylindrically  symmetric  heat  flow  equation  for  a  mixture  of  two 
energetic  materials  with  a  time  dependent  boundary  temperature.  Our  work  is  similar 
in  many  respects  to  recent  papers  by  Drake  [2],  and  McGuire  and  Tarver  [3],  which 
model  time  to  explosion  in  One  Dimensional  Time  to  Explosion  experiments  (ODTX), 
but  there  are  important  differences.  The  ODTX  experiments  impose  a  constant 
boundary  temperature  at  the  surface  of  a  heavily  confined,  spherical  explosive  sample, 
whereas  the  SSCB  test  applies  a  steadily  increasing  temperature  to  the  surface  of  a 
more  lightly  confined  explosive,  and  tests  are  typically  conducted  at  both  fast 
(approximately  l°C/second)  and  slow  (approximately  0.1°C/ second)  heating  rates. 

We  implemented  the  model  in  two  stages.  In  the  first  stage  we  used  a  zero  order 
kinetic  scheme  with  an  Arrhenius  temperature  dependent  rate  coefficient  to  describe 
the  decomposition  of  the  explosives,  and  we  assumed  that  the  specific  heat  and 
thermal  conductivity  of  the  compositions  were  independent  of  temperature.  This 
simple  model  was  sufficient  to  explain  the  trends  in  the  data  observed  by  Dagley  et  al., 
but  in  general  the  simulated  results  were  lower  than  the  experimental  values  for  both 
time  to  ignition  and  the  surface  temperature  at  ignition. 

We  then  considered  several  refinements  to  the  model,  including  a  three  term  kinetic 
scheme  for  the  decomposition  of  the  RDX,  and  allowance  for  the  temperature 
dependence  of  the  thermal  conductivity  and  specific  heat  of  the  materials.  Where 
appropriate,  these  were  included  in  the  calculations,  and  generally  lead  to  improved 
agreement  of  the  simulated  results  with  experiment. 

No  attempt  has  been  made  at  this  stage  to  model  the  violence  of  the  reaction  event. 
Current  understanding  of  the  physical  processes  involved  in  the  transition  to  either 
deflagration  or  detonation  in  a  cookoff  event  is  vague  [4],  although  encouraging 
progress  on  the  problem  has  been  made  recently  by  Cook  and  Haskins  [5]. 
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2.  Simplified  Model 


Dagley  et  al.  [1]  chose  RDX/Elvax  210  (95:5)  as  the  reference  composition  and  from 
this  mixtures  containing  TATB  and  PETN  at  levels  ranging  from  5%  to  35%  were 
prepared  (all  samples  were  pressed  to  90%  TMD),  and  their  cookoff  response  was 
determined  in  the  Super  Small-scale  Cookoff  Bomb.  The  SSCB  has  been  described  in 
detail  by  Parker  [6],  and  a  recent  report  by  Jones  and  Parker  [7]  described  numerical 
simulations  of  heat  flow  in  the  related  Small-scale  Cookoff  Bomb.  One  of  the 
conclusions  of  that  report  was  that  the  temperature  distribution  within  the  SCB  coidd 
be  accurately  described  using  a  one-dimensional  heat  flow  model,  and  so  we  have 
adopted  a  similar  one-dimensional  model  for  the  SSCB. 


2.1  Numerical  Solution 

The  radial  heat  flow  in  the  SSCB  for  a  multi-component  explosive  mixture  is  described 
by  the  following  equation  [8]: 


rdT  it 


jdvT 

{  dr 


1  dT 

~  H - 

2  r  dr 


(1) 


where  p  is  density,  C  is  specific  heat,  X  is  the  thermal  conductivity  of  the  mixture,  T  is 
temperature,  t  is  time,  and  r  is  radial  distance.  S  describes  the  rate  of  heat  generation 
per  unit  volume  at  temperature  T.  If  we  neglect  reactant  depletion  (ie.  zero  order 
kinetics)  and  assume  that  the  rate  of  reaction  of  each  of  the  species  varies  in 
accordance  with  the  Arrhenius  equation  then  S  has  the  form 


S  =  f>i"sQ.Ae*p[-£s/Kd  (2) 

S=1 

where  ms  is  the  mass,  Q.  the  heat  of  reaction  per  unit  mass,  As  the  pre-exponential 
factor,  Es  the  activation  energy,  and  NS  is  the  total  number  of  species.  R  is  the  gas 
constant. 


Equation  (1)  is  solved  using  an  operator  splitting  technique.  First  the  source  term  is 
set  to  zero  and  the  equation  is  discretised  using  a  standard  Forward  Time  Centred 
Space  (FTCS)  scheme  [9].  This  results  in  the  following  explicit  equation  for  the 
temperature  at  the  (n+l)th  time  step  at  node  i  in  terms  of  the  temperatures  at  the  nth 
time  step  at  nodes  i-1,  i,  and  r+1 


rpli  i  1  _  j'li  ,  C(.At 
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where  a  is  the  thermal  diffusivity  and  is  defined  by 


Equation  (3)  does  not  apply  on  the  axis  of  the  cylinder  where  z=l.  For  this  special  case 
we  apply  L'Hopital's  rule  to  the  third  term  on  the  left  side  of  equation  (1)  to  show  that 


Equation  (1)  then  becomes 


1  dT  d2T 

- > — tt  as  r— > 0 

r  dr  dr2 


dT  d2T 
pc*-_a“s 


and  has  a  finite  difference  representation  of  the  form 


7 f+i  =  T-1  +  {2 T?+1  ~  27}"} 


where  continuity  of  the  derivative  on  the  axis  has  also  been  used  in  deriving  the  finite 
difference  form.  Solutions  to  equations  (3)  and  (6)  are  stable  provided  that  the  time 
step  At  and  spatial  step  Ar  satisfy  the  following  inequality 

£M<I  (7) 

(Ar)2  2 

In  the  simulation  results  presented  here  At  was  always  chosen  to  be  considerably 
less  than  the  minimum  value  imposed  by  equation  (7),  and  several  values  of  At  and  Ar 
were  used  to  ensure  that  the  results  were  independent  of  the  spatial  and  time  step 
resolution. 

Solution  of  equation  (1)  proceeds  by  updating  the  temperature  at  each  node  point  at 
the  end  of  a  time  step  At  using  either  equation  (3)  (for  zVT)  or  equation  (6)  (for  z=l). 
The  special  case  z=N  corresponds  to  the  surface  of  the  SSCB.  The  temperature  at  this 
point  is  obtained  from  a  look-up  table  which  lists  the  experimentally  measured 
temperature  as  a  function  of  time  at  the  surface  of  the  SSCB  for  both  fast  and  slow 
cookoff.  These  values  were  determined  during  calibration  runs  using  an  SSCB  with  an 
inert  material  whose  thermal  properties  closely  matched  those  of  the  explosive  fills. 
The  temperature  increase  due  to  the  decomposition  of  the  explosive  is  then  included 
by  updating  the  temperature  at  each  node  point  from  the  source  term  S  using  the 
expression 


1SJ 

1 T1  -  Tr'+^msQsAsexP[-EjRTr\ 
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RDX  melts  at  a  temperature  of  204°C,  and  as  the  surface  temperature  at  ignition 
under  fast  cookoff  conditions  for  RDX/Elvax  210  (95:5)  is  240°C  an  appreciable 
fraction  of  the  RDX  content  in  the  SSCB  will  be  present  in  the  liquid  state.  Similar 
comments  apply  for  slow  cookoff,  where  the  surface  temperature  at  reaction  is  218°C. 
Hence  it  is  important  to  include  the  effect  of  melting.  We  did  this  using  the  following 
procedure;  as  soon  as  the  temperature  at  a  particular  node  reached  the  melting  point 
of  the  material  we  held  the  temperature  at  that  node  constant  at  the  melting  point 
value  until  an  energy  equivalent  to  the  latent  heat  of  fusion  of  the  material  had  been 
absorbed.  In  the  simple  model  presented  in  this  section  we  allowed  only  the  RDX 
component  to  undergo  melting.  PETN  however  has  a  lower  melting  point  of  140°C, 
and  we  have  allowed  for  melting  of  both  components  in  a  two  component  mixture  in  a 
later  section. 


2.2  Material  Parameters 


Solution  of  equation  (1)  for  a  two  component  explosive  mixture  requires  specification 
of  the  thermal  conductivity,  specific  heat,  density,  heat  of  reaction,  pre-exponential 
factor  and  activation  energy  for  each  of  the  component  species.  In  the  simplified 
model  used  here  we  have  made  the  approximation  that  the  value  of  each  of  these 
material  properties  is  independent  of  temperature.  Choosing  appropriate  values  for 
each  of  these  parameters  is  difficult  however  as  there  are  considerable  discrepancies 
between  some  of  the  values  quoted  in  the  literature. 

The  thermal  conductivity  (at  room  temperature)  of  RDX  is  given  by  Rogers  [10]  as 
2.53  x  10"4  cal  cm4  s'1  °C~1,  by  Zinn  and  Mader  [11]  as  7  x  10'4  cal  cnr1  s'1  “C"1,  and 
by  McGuire  and  Tarver  [3]  and  Drake  [2]  as  6.22  x  10'4  cal  cm"1  s"1  °C_1.  There  is  less 
disagreement  regarding  the  value  of  the  specific  heat  of  RDX  at  room  temperature; 
Dobratz  and  Crawford  [12]  quote  a  value  of  0.26  x  10"4  cal  g"1  °C'1,  while  McGuire  and 
Tarver  use  0.24  x  10"4  cal  g"1  °C"1.  There  are  also  discrepancies  in  the  values  for  the 
activation  energy,  pre-exponential  factor,  and  heat  of  reaction  for  RDX.  Zinn  and 
Mader  use  E  =  47.5  kcal  M"1,  A  =  3.162  x  1018  s"1  and  Q  =  500  cal  g"1,  Rogers  uses  E  = 
47.1  kcal  M'1  ,  A  =  2.02  x  1018  s'1  and  Q  =  500  cal  g'1,  while  Hutchinson  [13]  uses  E  = 
48.5  kcal  M"1,  A  =  8.6  x  1018  s"1  and  Q  =  549  cal  g"1. 

The  values  we  have  chosen  for  RDX  are  shown  in  Table  1.  These  were  determined  by 
performing  a  number  of  numerical  simulations  of  the  time  to  reaction  and  temperature 
at  reaction  for  RDX  for  both  fast  and  slow  cookoff  events,  and  then  chosing  those 
values  which  gave  the  best  overall  fit  to  the  experimental  results.  The  values  chosen 
for  specific  heat  and  thermal  conductivity  also  have  the  advantage  of  agreeing  with 
the  temperature  dependent  values  of  these  quantities  at  room  temperature,  which  will 
facilitate  comparisons  with  experiment  in  a  later  section  where  we  specifically  include 
the  effect  of  the  temperature  dependence  of  these  properties. 
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Table  1:  Thermochemical  Constants 


RDX 

PETN 

TATB 

Density  (g.ml  "^) 

1.80 

1.74 

1.84 

Thermal  conductivity  (cal  cm'k'^  °C“^) 

0.00062 

0.0006 

0.00191 

Specific  Heat  (cal  g"^°C'^) 

0.24 

0.272 

0.26 

> 

1 

t-1 

3.162X1018 

6.300X1019 

3.180xl019 

Q  (cal  g'1) 

500 

300 

600 

E  (kcal  M"1) 

48.0 

47.0 

59.9 

Latent  Heat  of  Fusion  (cal  g'^) 

38.4 

10.0 

- 

Melting  Point  (°C) 

204 

141 

320 

The  values  used  for  PETN  and  TATB  have  been  taken  primarily  from  the  paper  by 
Rogers  [10],  although  the  specific  heat  and  thermal  conductivity  values  have  been 
taken  from  the  more  recent  papers  by  McGuire  and  Tarver  [3],  and  Drake  [2]. 
Calculations  were  also  performed  at  both  crystal  density  and  90%  TMD,  and  it  was 
found  that  the  slight  change  in  density  had  a  negligible  effect  on  the  results. 

Table  2  shows  the  calculated  time  to  explosion  and  surface  temperature  at  reaction 
for  RDX/Elvax  210  (95:5)  for  both  fast  and  slow  cookoff  conditions  using  the  values 
for  the  physical  constants  shown  in  Table  1  (the  values  in  brackets  are  the  percentage 
deviation  from  the  experimental  results).  The  effect  of  the  Elvax  210  on  the  behaviour 
of  RDX  has  been  neglected  at  this  stage,  and  we  will  discuss  this  point  in  more  detail 
in  a  later  section.  Thermal  explosion  was  considered  to  have  occured  when  the 
temperature  at  some  point  in  the  explosive  reached  a  predetermined  high  value.  For 
the  calculations  reported  here  this  value  was  550  °C,  although  the  actual  value  used 
had  negligible  effect  on  the  results.  The  calculations  were  performed  both  with  and 
without  the  inclusion  of  melting,  and  the  results  clearly  show  the  importance  of 
including  this  effect,  particularly  for  the  slow  cookoff  results. 


Table  2:  Simulated  results  for  RDX/Elvax  210  (95:5) 


FAST  COOKOFF 

Melting  Neglected 

Melting  Included 

Experiment 

Time  to  Event  (s) 
Surface  Temp.  °C 

242  (1.6%) 

232  (3.3%) 

256  (3.6%) 

241  (0.4%) 

246 

240 

SLOW  COOKOFF 

Melting  Neglected 

Melting  Included 

Experiment 

Time  to  Event  (s) 
Surface  Temp.  °C 

1208  (27%) 

195  (11%) 

1461  (11.7%) 

210  (3.6%) 

1654 

218 
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The  results  shown  in  Table  2  will  be  improved  in  a  later  section  when  we  consider  a 
more  accurate  model  for  RDX,  but  they  are  adequate  as  a  starting  point  for  modelling 
the  effect  of  the  addition  of  varying  amounts  of  PETN  and  TATB  on  the  cookoff 
behaviour  of  RDX/Elvax  210. 

To  calculate  the  thermal  conductivity  of  the  explosive  mixtures  we  have  used  an 
expression  originally  derived  by  Maxwell  [14],  and  more  recently  by  Jeffrey  [15],  who 
derived  an  expression  for  the  thermal  conductivity  to  a  higher  order  in  the  species 
concentration.  Jeffrey  studied  the  conduction  of  heat  through  a  stationary,  random 
and  statistically  homogeneous  suspension  of  spherical  particles  in  a  matrix  of  uniform 
conductivity  under  the  condition  that  the  volume  fraction  of  the  particles  is  small.  If 
the  matrix  has  thermal  conductivity  At,  and  the  spheres  a  thermal  conductivity  A,, 
then  to  first  order  in  the  concentration  c  the  effective  conductivity  A  is  given  by 

A  =  A1{l  +  3[k}  (9) 

where 

Ma2-x,)/(a2  +  x,)  (io) 

Equation  (9)  is  valid  when  the  concentration  is  small  enough  to  make  all  interactions 
between  the  spheres  negligible.  Jeffrey  has  extended  the  expression  for  A  to  second 
order  in  c  by  allowing  for  interactions  between  pairs  of  spheres,  but  in  keeping  with 
the  simplicity  of  the  model  adopted  for  each  of  the  explosive  components  we  have 
calculated  the  thermal  conductivity  of  the  mixture  using  the  lower  order  correction 
only.  Use  of  equation  (9)  gives  results  which  are  very  close  to  those  calculated  by 
assuming  that  the  thermal  conductivity  of  the  mixture  is  simply  given  by  the  mass 
average  of  the  values  for  the  individual  explosives,  which  is  the  method  used  by 
McGuire  and  Tarver  in  their  calculations  [3]. 


3.  Results  of  Simplified  Model 


As  described  by  Dagley  et  al  [1],  TATB  and  PETN  were  added  to  the  RDX/Elvax  210 
(95:5)  reference  composition  to  modify  its  cookoff  response  by  incorporating 
explosives  having  different  thermal  stabilities.  TATB  has  a  much  higher  critical 
temperature  than  RDX  and  was  expected  to  act  as  an  inert  material  during  the  heating 
phase,  thus  lowering  the  explosive  output  during  accidental  thermal  stimulus,  but 
retaining  full  explosive  output  when  initiated  in  the  normal  maimer.  PETN  has  a 
lower  critical  temperature  and  was  expected  to  lead  to  reaction  at  lower  temperatures 
and  to  produce  milder  responses  due  to  early  release  of  confinement  resulting  from 
the  initial  reaction  of  the  PETN. 

At  the  slow  heating  rate  both  the  PETN  and  the  TATB  succeeded  in  moderating  the 
violence  of  the  cookoff  reaction.  As  the  PETN  content  increased  the  explosive  surface 
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temperature  at  reaction  decreased  and  milder  responses  were  also  obtained  in  most 
cases.  Figure  1  shows  the  experimental  results  and  the  predictions  of  the  model  for  the 
surface  temperature  at  reaction  as  a  function  of  PETN  content  for  the  slow  heating 
rate.  (In  all  the  Figures  presented  here  the  open  circles  represent  the  experimental 
results  and  the  filled  circles  represent  the  simulated  results.  The  dashed  line  through 
the  experimental  data  typically  denotes  a  second  order  polynomial  fit  to  the  data,  and 
the  solid  line  a  second  order  polynomial  fit  to  the  model  results.)  As  noted  in  the 
previous  section,  the  simple  model  employed  here  uses  temperature  independent 
values  for  the  specific  heat  and  thermal  conductivity  for  each  explosive,  and  the 
particular  values  chosen  give  a  better  fit  to  the  data  for  RDX/Elvax  210  (95.5)  at  the 
fast  heating  rate  rather  than  at  the  slow  heating  rate.  This  is  the  cause  of  the  offset 
between  the  two  curves  in  Figure  1.  Apart  from  this  discrepancy,  the  model  provides  a 
good  qualitative  explanation  of  the  trend  in  the  experimental  data.  As  expected,  the 
decrease  in  surface  temperature  at  reaction  is  due  to  the  PETN  content  undergoing 
thermal  runaway  before  the  RDX  begins  to  react  appreciably.  If  the  RDX  source  term 
is  omitted  from  equation  (2)  then  almost  identical  results  are  obtained,  indicating  that 
the  thermal  runaway  is  triggered  by  the  PETN.  Figure  2  shows  the  time  to  reaction  as 
a  function  of  PETN  content  for  the  slow  heating  rate  and  similar  results  are  obtained. 
Again  the  decrease  in  time  to  reaction  as  the  percentage  of  PETN  in  the  composition  is 
increased  is  qualitatively  reproduced  by  the  model. 

For  different  levels  of  TATB  at  the  slow  heating  rate  the  reactions  generally  occurred 
in  the  same  temperature  range  as  the  RDX/Elvax  210,  indicating  that  the  reaction  was 
being  triggered  by  the  RDX,  and  the  cookoff  responses  were  generally  milder  than  for 
the  reference  composition.  Figure  3  shows  the  experimental  results  and  the  predictions 
of  the  model  for  the  surface  temperature  at  reaction  as  a  function  of  TATB  content  for 
the  slow  heating  rate.  The  experimental  results  show  a  very  slight  increase  in  surface 
temperature  at  reaction  as  the  TATB  content  is  increased  and  the  model  accurately 
reproduces  this  trend,  with  the  rate  of  increase  in  temperature  being  very  similar  for 
both  curves.  Again  there  is  an  approximately  constant  offset  between  the  two  lmes 
because  of  the  nature  of  the  model  used  for  the  RDX/Elvax  210  (95.5)  composition. 
The  very  small  effect  which  the  TATB  has  on  the  surface  temperature  at  reaction  is 
due  to  the  very  slight  change  in  thermal  conductivity  of  the  mixture.  The  TATB  is 
thermally  very  stable  in  the  temperature  range  shown  in  Figure  3  and  the  reaction  is 
triggered  by  the  RDX.  Removal  of  the  TATB  source  term  in  equation  (2)  leads  to 
exactly  the  same  results. 

Figure  4  shows  the  time  to  reaction  as  a  function  of  TATB  content  for  the  slow 
heating  rate.  The  model  predicts  a  steady  increase  in  time  to  reaction  as  the  TATB 
content  increases  whereas  a  straight  line  fit  to  the  experimental  data  indicates  a 
decrease  in  time  to  reaction  as  TATB  content  increases.  There  is  considerable  scatter  m 
the  experimental  data  shown  in  Figure  4  however  and  more  data  is  required  to  get  a 
true  representation  of  the  experimental  trend. 
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Figure  1:  Surface  temperature  at  reaction  as  a  function  of  PETN  content.  The  open  circles 
are  the  experimental  points  and  the  filled  circles  are  the  predictions  of  the  simple  model. 
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Figure  2:  Time  to  reaction  as  a  function  of  PETN  content.  The  open  circles  are  the 
experimental  points  and  the  filled  circles  are  the  predictions  of  the  simple  model. 


DSTO-TR-0090 


SSCB  Tests  of  RDX/TATB/EIvax  210 

SLOW  HEATING  RATE 


%TATB 


Figure  3:  Surface  temperature  at  reaction  as  a  function  ofTATB  content.  The  open  circles 
are  the  experimental  points  and  the  filled  circles  are  the  predictions  of  the  simple  model. 


Figure  4:  Time  to  reaction  as  a  function  of  TATB  content.  The  open  circles  are  the 
experimental  points  and  the  filled  circles  are  the  predictions  of  the  simple  model. 


9 


DSTO-TR-0090 


At  the  fast  heating  rate  addition  of  PETN  had  no  appreciable  effect  on  the  reaction 
temperature  until  a  considerable  amount  of  PETN  (25%-35%)  had  been  added,  and  all 
the  compositions  gave  more  violent  responses  than  the  reference  composition. 
Figure  5  shows  the  experimental  results  and  the  predictions  of  the  model  for  the 
surface  temperature  at  reaction  as  a  function  of  PETN  content  for  the  fast  heating  rate. 
The  experimental  results  are  again  quite  scattered,  but  the  overall  trend  is  to  a  slow 
decrease  in  surface  temperature  at  reaction  for  small  percentages  of  PETN,  followed 
by  a  more  rapid  decrease  at  higher  concentrations.  The  model  however  predicts  the 
reverse,  with  the  maximum  decrease  in  surface  temperature  occurring  at  low 
concentrations  of  PETN  and  then  the  temperature  appears  almost  to  level  off  for 
higher  concentrations.  Similar  results  are  shown  in  Figure  6,  where  we  plot  the  time  to 
reaction  as  a  function  of  PETN  content  at  the  fast  heating  rate.  At  low  PETN 
concentrations  the  two  curves  show  opposing  trends,  with  the  experimental  values 
increasing  with  concentration  and  the  model  predictions  decreasing  with 
concentration.  These  conflicting  results  may  be  due  to  some  of  the  assumptions  made 
in  the  simple  model,  particularly  the  assumption  of  temperature  independent  values 
for  the  thermal  conductivity  and  specific  heat.  We  will  pursue  this  point  further  in  the 
next  section. 

The  addition  of  TATB  had  little  effect  on  either  reaction  temperature  or  response  at 
the  fast  heating  rate;  reaction  continued  to  occur  in  the  same  temperature  range  as  for 
the  RDX/Elvax  210,  and  the  responses  were  generally  similar  to  those  of  the 
RDX/Elvax  210.  Figure  7  shows  the  experimental  results  and  the  predictions  of  the 
model  for  the  surface  temperature  at  reaction  as  a  function  of  TATB  content  for  the 
fast  heating  rate.  The  experimental  results  are  fairly  scattered  but  show  a  trend 
towards  a  very  slight  increase  in  surface  temperature  at  reaction  as  the  TATB  content 
is  increased.  The  model  similarly  predicts  an  increase  in  temperature  with  increasing 
TATB  content. 

Figure  8  shows  the  time  to  reaction  as  a  function  of  TATB  content  for  the  fast  heating 
rate.  The  experimental  data  is  extremely  scattered,  but  shows  a  trend  towards 
increasing  temperatures  with  increasing  TATB  content.  This  is  also  predicted  by  the 
model,  although  the  rate  of  increase  is  much  slower  than  that  shown  experimentally. 
The  discrepancy  between  the  results  is  probably  a  combination  of  insufficient 
experimental  data  and  the  inadequacy  of  some  of  the  assumptions  employed  in  the 
model. 
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SSCB  Tests  of  RDX/PETN/Elvax  210 


Figure  5:  Surface  temperature  at  reaction  as  a  function  of  PETN  content.  The  open  circles 
are  the  experimental  points  and  the  filled  circles  are  the  predictions  of  the  simple  model. 


SSCB  Tests  of  RDX/PETN/Elvax  210 


Figure  6:  Time  to  reaction  as  a  function  of  PETN  content.  The  open  circles  are  the 
experimental  points  and  the  filled  circles  are  the  predictions  of  the  simple  model. 
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Figure  7:  Surface  temperature  at  reaction  as  a  function  ofTATB  content.  The  open  circles 
are  the  experimental  points  and  the  filled  circles  are  the  predictions  of  the  simple  model. 
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Figure  8:  Time  to  reaction  as  a  function  of  TATB  content.  The  open  circles  are  the 
experimental  points  and  the  filled  circles  are  the  predictions  of  the  simple  model. 
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4.  Improved  Model 


Since  the  early  work  of  Rogers  [10],  Zinn  and  Mader  [11],  and  Zinn  and  Rogers  [16], 
more  refined  models  of  the  thermal  behaviour  of  many  explosives  have  appeared  in 
the  literature.  Prominent  among  these  is  the  work  of  Catelano  et  al.  [17],  Tarver  et  al. 
[18],  and  McGuire  and  Tarver  [3].  These  authors  have  considered  thermal 
decomposition  of  confined  explosives  in  one-dimensional  geometries  and  have 
obtained  good  agreement  between  experimental  results  and  calculations  in  One- 
Dimensional-Time-to-Explosion  (ODTX)  tests.  The  calculated  results  include  the 
temperature  dependence  of  the  specific  heat  and  thermal  conductivity,  as  well  as 
multi- term  schemes  for  the  decomposition  reactions.  In  this  section  we  implement  the 
model  for  RDX  described  by  McGuire  and  Tarver  [3]  and  consider  the  effect  of  these 
improvements  on  the  calculated  time  to  explosion  and  surface  temperature  for  the 
SSCB  test  on  RDX/Elvax  210. 


4.1  Temperature  Dependant  Thermal  Properties 


When  the  thermal  conductivity  and  specific  heat  are  temperature  dependant,  the 
radial  heat  flow  in  the  SSCB  is  described  by  the  following  equation  [19] 


dt  dr  l  dr  J  r  dr 


The  finite  difference  form  of  equation  (11)  is  found  by  replacing  each  of  the  differential 
terms  by  the  following  expressions  [20] 


dt 


3/  ^ -T  ■  \  fl  (  rr>J1  rj^Tl  \  f  rjift  rpfl? 

fl/T\C'J  I  .VI  1M~1i  VI  h  -fl-1 

*{m *  J,  J ‘ v 


2 X{T)  dTY  2X1  tu\ -1 

r  dr ).  (z-l)Ar  2A r 


^  rr,ii 

!alV]  ...  for 


-2V] . for  i-1 

y  dr  Avz  >-  J 


where  Xnhin  and  X']  V2  are  defined  by 
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^1l/2=i(X','>+X’,)  'and  -  ^  V2=  +  l).  (16) 

which  results  in  the  following  expression  for  the  temperature  at  the  (n+l)th  time  step 
at  node i 


for  i*l,  and 


L-^-\T"x-T"} . for  i=l 

"  (aA2  I  1  +  1  1  -I 


pc’’  (A r) 


Solution  of  equation  (11)  proceeds  in  an  iterative  fashion.  We  begin  with  an  initial 
value  for  the  temperature  at  all  grid  points.  Values  for  the  thermal  conductivity  and 
specific  heat  are  then  calculated  and  equations  (17)  and  (18)  solved  to  get  the  new 
temperature  values.  New  values  for  the  thermal  conductivity  and  specific  heat  are 
then  calculated  using  the  new  temperature  values  and  then  equations  (17)  and  (18) 
are  solved  again  to  get  new  estimates  for  the  temperature.  The  cycle  is  repeated  until 
successive  changes  in  temperature  fall  below  some  prescribed  small  limit.  The 
temperature  increase  from  the  source  term  is  then  added  in  the  same  manner  as 
before,  and  the  cycle  is  then  repeated  for  the  next  time  step. 


We  have  also  included  the  effect  of  reactant  depletion  at  this  stage  by  writing  the 
source  term  as 


S  =  p(Moi>exp[-E/RT] 


(19) 


where  0)  is  the  mass  fraction  of  undecomposed  RDX.  The  time  dependence  of  co  is 
described  by 

^  =  -Acoexp[-£/RT]  (20) 

ot 

As  an  illustration  of  the  effect  of  the  inclusion  of  temperature  dependent  thermal 
properties  and  reactant  depletion,  the  time  to  explosion  and  surface  temperature  have 
been  calculated  for  RDX  in  the  SSCB  for  both  fast  and  slow  cookoff.  The  temperature 
dependence  of  the  specific  heat  and  thermal  conductivity  were  taken  from  the  paper 
by  McGuire  and  Tarver  [3].  Over  the  temperature  range  considered  the  thermal 
conductivity  of  the  RDX  decreases  by  approximately  30%  to  40%  as  the  temperature 
increases,  while  the  specific  heat  increases  by  approximately  50%.  The  results  are 
presented  in  Table  3  (the  experimental  values  are  for  the  reference  composition, 
RDX/Elvax  210  (95:5),  not  pure  RDX).  These  figures  show  that  the  inclusion  of  both 
these  effects  has  a  marked  effect  on  the  calculated  times  and  temperatures  for  the 
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SSCB  cookoff  tests.  Agreement  with  experiment  has  been  improved  considerably  for 
slow  cookoff,  and  even  the  fast  cookoff  results  have  improved  slightly. 


Table  3:  Simulated  results  for  RDX/Elvax  210  (95:5)  using  Improved  Model 


FAST 

COOKOFF 

X  and  C  fixed 

X  and  C  variable 

Experiment 

Time  to  Event  (s) 

258  (5%) 

256  (4%) 

246 

Surface  Temp.  °C 

242  (1%) 

241  (1%) 

240 

SLOW 

COOKOFF 

X  and  C  fixed 

X  and  C  variable 

Experiment 

Time  to  Event  (s) 

1400  (15%) 

1645  (0.5%) 

1654 

Surface  Temp.  °C 

206  (6%) 

216  (1%) 

218 

Figures  9  and  10  show  the  radial  temperature  profile  within  the  SSCB  at  selected 
times  for  both  fast  and  slow  cookoff.  Under  fast  cookoff  conditions  (Figure  9)  there  is  a 
considerable  thermal  gradient  within  the  SSCB,  the  temperature  is  highest  at  the  edge 
of  the  cookoff  bomb,  and  thermal  runaway  commences  very  close  to  the  surface  of  the 
explosive.  For  slow  cookoff  (Figure  10)  there  is  a  more  uniform  temperature 
distribution  within  the  SSCB  and  the  reaction  commences  closer  to  the  centre  of  the 
bomb.  Figure  10  also  indicates  that  an  appreciable  fraction  of  the  RDX  is  present  in 
the  liquid  state  for  approximately  two  minutes  prior  to  the  start  of  reaction.  Recent 
work  using  a  more  extensive  thermocouple  array  within  the  SSCB  has  confirmed  that 
melting  of  the  RDX  does  occur  at  the  centre  of  the  SSCB  approximately  100  seconds 
prior  to  reaction,  although  the  extent  of  the  melting  is  less  than  that  shown  in  Figure 
10  [25]. 


4.2  Multi-term  Kinetic  Scheme 

The  chemical  decomposition  model  for  RDX  described  by  McGuire  and  Tarver  [3] 
consists  of  three  reactions,  represented  by  the  following  sequence: 

A  — >  B  — >  2C  — »  D  (21) 

where  A  represents  RDX,  B  represents  H2C=N-N02,  C  represents  (CH2+N2)  or 
(HCN+HN02),  and  D  represents  the  final  products.  As  the  decomposition  proceeds 
the  concentration  of  each  of  the  species  is  described  by  the  following  set  of  coupled 
rate  equations: 
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dA  ,  A 

—  =  -LA 
dt  1 

(22) 

—  -  LA  -  k2B 

(23) 

dt 

~  =  k2B-  k3C2 
dt  2  3 

(24) 

—  =  t3C2 

(25) 

dt 

where  the  rate  constants  kj  are  defined  by  kj=Zjexp(-Ej/RT).  The  source  term  in  the 
heat  flow  equation  for  this  model  has  the  form 

NAq]Z,e  W  +  NBq2Z2e  WRT  +  N2cq3Z3e^RT  (26) 

where  and  N|j  are  the  respective  mole  fractions  defined  such  that 

Na+Nj3+Nc+N|)=1. 

RDX/Elvax  210 


Figure  9 :  Simulated  radial  temperature  profile  within  the  SSCB  at  selected  times . 
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Figure  10:  Simulated  radial  temperature  profiles  within  the  SSCB  at  selected  times. 


The  above  kinetic  scheme  was  included  in  the  computer  code  described  in  section 
4.1,  which  was  then  used  to  model  both  fast  and  slow  cookoff  for  RDX/Elvax  210  in 
the  SSCB.  The  results  are  shown  in  Table  4.  These  indicate  that  the  more  sophisticated 
decomposition  scheme  provides  worse  agreement  with  experiment  than  does  the 
simple  first  order  scheme.  The  reason  for  this  is  believed  to  be  related  to  the  differing 
amounts  of  confinement  in  the  SSCB  and  ODTX  tests.  The  papers  by  Tarver  et  al.  [18] 
and  Catalano  et  al.  [17]  have  shown  that  the  time  to  explosion  is  strongly  influenced 
by  the  void  volume  of  the  containment  vessel.  Even  the  addition  of  a  relatively  small 
void  space  in  the  ODTX  experiments  resulted  in  almost  a  doubling  of  the  time  to 
explosion  for  LX-04  (85%  HMX,  15%  binder).  In  the  earlier  experimental  work  of 
Rogers  [10],  Zinn  and  Mader  [11],  and  Zinn  and  Rogers  [16],  the  explosive  samples 
were  less  well  confined  and  the  explosive  decomposition  products  were  easily  vented, 
whereas  in  the  ODTX  experiments  great  care  was  taken  to  ensure  that  none  of  the 
gaseous  products  were  allowed  to  escape  from  contact  with  the  high  explosive 
sample.  Calculations  of  critical  temperatures  for  the  explosives  TNT,  TATB  and  LX-10 
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using  a  multi-term  decomposition  scheme  similar  to  the  one  above  gave  good 
agreement  with  results  from  ODTX  tests,  but  were  appreciably  lower  than  values 
calculated  using  the  earlier  single  step  schemes  [18].  Similar  reductions  in  surface 
temperature  at  reaction  are  observed  here  for  the  RDX/Elvax  210  composition  shown 
in  Table  4.  The  SSCB  contains  a  relatively  large  void  space  at  the  top  of  the  cylinder, 
and  the  top  confining  plate  contains  a  small  hole  which  allows  access  for  the 
thermocouple,  and  remains  unsealed  during  the  test.  Hence  the  gaseous 
decomposition  products  are  relatively  easily  vented,  and  the  SSCB  is  more  aptly 
described  by  the  single  step  reaction  schemes  constructed  to  model  the  earlier 
experiments.  When  the  multi-term  kinetic  scheme  is  used,  as  in  the  results  shown  in 
Table  4,  the  effect  is  to  reduce  both  the  time  to  explosion  and  the  surface  temperature 
at  explosion,  for  both  fast  and  slow  cookoff.  A  modified  form  of  SSCB  test  is  to  be 
constructed  shortly  at  MRL.  This  will  include  pressure  transducer  diagnostics,  and  the 
cylinder  will  be  sealed  to  ensure  that  gaseous  decomposition  products  will  be 
contained  until  the  cylinder  is  ruptured  [21].  Under  these  conditions  the  SSCB  results 
are  expected  to  agree  more  closely  with  the  calculated  values  shown  in  Table  4. 


Table  4:  Simulated  results  for  RDX/Elvax  210  (95:5)  using  Improved  Model 


FAST 

COOKOFF' 

First  order  kinetics 

Multi  term  kinetics 

Experiment 

Time  to  Event  (s) 

263  (7%) 

235  (4%) 

246 

Surface  Temp.  °C 

245  (2%) 

226  (6%) 

240 

SLOW 

COOKOFF 

First  order  kinetics 

Multi  term  kinetics 

Experiment 

Time  to  Event  (s) 

1580  (5%) 

1315(20%) 

1654 

Surface  Temp.  °C 

214  (2%) 

201  (8%) 

218 

4.3  Effect  of  Binder 


The  calculations  to  date  have  used  a  reference  composition  consisting  of  100%  RDX 
instead  of  the  actual  composition,  consisting  of  95%  RDX  and  5%  Elvax  210.  We  now 
consider  the  effect  of  this  approximation  on  the  accuracy  of  the  calculated  results.  Data 
on  the  thermophysical  constants  of  Elvax  210  are  difficult  to  obtain.  A  recent  paper  by 
Dagley  et  al.  [22]  quotes  the  thermal  conductivity  of  Elvax  265  (a  related  polymer 
which  should  have  very  similar  thermophysical  properties  to  Elvax  210)  as  7.14  x  10"4 
cal  cnrV1  0C_1  at  30°C  and  5.24  x  10'4  cal  arris'1  XT1  at  200°C.  These  are  very  close 
to  the  thermal  conductivity  values  for  RDX  used  by  McGuire  et  al.  [3]  of  6.22  x  10'4  cal 
cm4s4  °C‘l  at  20°C  and  4.85  x  10~4  cal  crn'ls'l  °C"1  at  160°C.  An  estimate  of  the 
effective  thermal  conductivity  of  the  RDX/Elvax  210  (95:5)  mixture  can  be  made  by 
assuming  that  the  composition  consists  of  spheres  of  RDX  uniformly  coated  with 
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Elvax  210.  Helsing  and  Grimvall  [23]  quote  the  following  expression  for  the  effective 
thermal  conductivity  of  such  a  composite  system  : 


where  phase  p  represents  the  RDX,  phase  a  represents  the  Elvax  210,  and  fa  and  /p 
represent  the  respective  volume  fractions.  When  this  expression  is  used  to  calculate 
the  thermal  conductivity  at  the  lower  temperature  a  value  of  6.10xl0'4  cal  cm'V1  °C~1 
is  obtained;  this  is  very  close  to  the  value  for  RDX  at  this  temperature,  indicating  that 
our  use  of  the  thermal  conductivity  of  RDX  for  the  effective  thermal  conductivity  of 
the  composite  system  is  quite  a  good  approximation. 

The  specific  heat  of  Elvax  210  is  quoted  by  Dagley  et  al.  [22]  as  being  0.57  cal 
g‘loC'1  at  150°C  and  0.63  cal  g'loC_1  at  220°C,  which  is  considerably  higher  than  the 
values  used  by  McGuire  et  al.  [3]  of  0.24  cal  g-l°C‘l  at  20°C  and  0.42  cal  g'4°C'^  at 
350°C  for  RDX.  The  effective  specific  heat  of  the  composite  material  will  still  be  very 
close  to  the  RDX  value  however  because  the  contribution  from  the  Elvax  210  (which  is 
directly  proportional  to  its  mass  fraction)  will  be  small  compared  with  the  contribution 
from  the  RDX.  Another  way  to  look  at  this  is  to  consider  that  the  higher  specific  heat 
of  the  Elvax  210  will  result  in  a  lower  thermal  diffusivity  for  this  material.  If  the 
thermal  diffusivity  is  calculated  using  X  from  equation  (27)  then  the  calculated  value  is 
very  close  to  the  thermal  diffusivity  of  RDX. 

The  discussion  above  is  based  on  the  use  of  partly  known  and  partly  inferred 
thermophysical  constants  for  the  binder.  Given  that  these  approximate  expressions 
result  in  estimates  for  the  effective  thermal  conductivity  and  specific  heat  for  the 
composite  material  which  are  very  close  to  those  for  RDX,  we  feel  justified  in 
continuing  to  omit  the  binder  from  the  calculations,  and  to  consider  the  reference 
composition  to  consist  entirely  of  RDX.  For  reasons  similar  to  those  discussed  above 
we  have  also  not  included  in  the  calculations  a  term  to  account  for  the  melting  of  the 
binder.  For  Elvax  210  this  occurs  around  60°C. 

While  we  expect  that  neglecting  the  thermophysical  properties  of  the  binder  will 
have  little  effect  on  the  thermophysical  properties  of  the  reference  composition,  and 
therefore  on  the  predicted  cookoff  temperature  and  time  to  event,  the  pressence  of  the 
binder  is  known  to  have  a  significant  effect  on  the  violence  of  the  reaction  event  [24]. 


5.  Results  of  Improved  Model  for  Explosive  Mixtures 


The  results  presented  in  Section  3  have  shown  that  a  simple  model  of  the  explosive 
mixtures  based  on  the  use  of  zero  order  kinetics  to  describe  the  reaction 
decomposition,  an  Arrhenius  temperature  dependence  for  the  rate  of  reaction. 
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temperature  independent  thermophysical  constants,  and  allowance  for  the  melting  of 
the  RDX,  is  capable  of  reproducing  the  observed  experimental  trends  in  the  time  to 
reaction  and  surface  temperature  at  reaction.  In  Section  4  we  showed  that  by  including 
reactant  depletion  (ie.  first  order  decomposition  reaction)  and  the  temperature 
dependence  of  the  specific  heat  and  thermal  conductivity  of  the  explosive  we  were 
able  to  obtain  improved  agreement  with  the  experimental  results  for  RDX/Elvax  210 
(95:5)  at  both  fast  and  slow  cookoff.  We  also  showed  that  the  decomposition  reaction 
schemes  devised  by  McGuire  and  Tarver  [3],  and  Drake  [2],  for  ODTX  experiments 
were  inappropriate  for  the  SSCB.  In  this  section  we  improve  our  models  for 
RDX/PETN/Elvax  210  and  RDX/TATB/Elvax  210  by  including  the  relevant 
refinements  described  in  Section  4. 

Our  models  for  the  two  explosive  mixtures  are  necessarily  different.  Data  on  the 
temperature  dependence  of  the  thermal  conductivity  and  specific  heat  of  PETN  were 
unavailable,  and  so  could  not  be  included.  PETN  however  melts  at  a  temperature  of 
approximately  141  °C  and  has  a  latent  heat  of  fusion  of  approximately  10  kcal/mole, 
and  so  allowance  for  the  melting  of  both  the  RDX  and  PETN  was  included  in  the 
model,  together  with  a  first  order  reaction  scheme  for  both  the  RDX  and  PETN.  The 
temperature  dependence  of  the  thermal  conductivity  and  specific  heat  of  TATB  has 
been  reported  by  both  McGuire  and  Tarver  [3]  and  Drake  [2],  and  was  included  in  the 
refined  model  for  RDX/TATB/Elvax  210.  First  order  reaction  schemes  were  employed 
for  both  RDX  and  TATB,  although  the  TATB  does  not  react  in  the  temperature  range 
considered  here.  No  allowance  for  TATB  melting  was  necessary,  since  its  melting 
point  (320°C  [12])  is  above  the  temperature  at  which  reaction  occurs  in  any  of  our 
experiments  or  calculations. 

Figure  11  shows  the  experimental  and  calculated  surface  temperature  at  reaction  as  a 
function  of  PETN  content  at  the  slow  heating  rate.  Considerable  improvement  over  the 
results  from  the  simple  model  (  Figure  1)  are  observed;  there  is  now  good  agreement 
with  experiment  for  an  appreciable  range  of  PETN  content,  although  the  calculated 
and  experimental  results  still  diverge  at  30%  to  40%  PETN  content.  This  may  be  due  to 
the  neglect  of  the  temperature  dependence  of  the  specific  heat  and  thermal 
conductivity  of  the  PETN,  or  simply  the  large  scatter  in  the  experimental  points.  If  the 
high  temperature  result  at  35%  PETN  content  was  in  error  and  the  true  value  was 
closer  to  the  low  temperature  result  then  the  calculated  and  experimental  curves 
would  be  in  agreement  over  the  complete  range.  Figure  12  shows  that  a  similar  degree 
of  improvement  can  be  seen  in  the  time  to  reaction  as  a  function  of  PETN  content  at 
the  slow  heating  rate. 

At  the  fast  heating  rate  the  improved  model  for  RDX/PETN/Elvax  210  shows 
considerably  better  agreement  with  experiment  than  the  simple  model.  Figure  13 
shows  the  surface  temperature  at  reaction  as  a  fimction  of  PETN  content.  The  simple 
model  (Figure  5)  predicted  a  rapid  decrease  in  temperature  for  small  amoimts  of 
PETN  followed  by  a  levelling  off  at  around  30%  to  40%  PETN  content,  while 
experiment  showed  that  the  temperature  decreased  only  slowly  at  small  percentages 
of  PETN,  and  then  dropped  quickly  in  the  20%  to  40%  range.  The  improved  model  has 
corrected  these  opposing  trends;  the  drop  in  temperature  is  now  quite  small  at  the 
lower  PETN  levels,  and  the  experimental  and  simulated  results  agree  well  to  about 
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20%  PETN  concentration.  For  higher  levels  of  PETN  content  the  simulated  result  is 
now  too  high,  but  this  may  again  be  the  result  of  the  neglect  of  the  temperature 
dependence  of  the  thermophysical  constants  for  PETN.  Even  better  agreement  can  be 
seen  in  Figure  14,  which  shows  the  time  to  reaction  as  a  function  of  PETN  content  at 
the  fast  heating  rate.  The  simple  model  (Figure  6)  showed  general  agreement,  but 
again  the  maximum  simulated  rate  of  decrease  in  temperature  occured  at  the  wrong 
end  of  the  range.  Figure  14  shows  excellent  agreement  between  the  simulated  curve 
and  a  straight  line  fit  to  the  experimental  data. 

For  RDX/TATB/Elvax  210  there  is  also  generally  better  agreement  with  experiment 
for  the  improved  model.  Figure  15  shows  surface  temperature  at  reaction  as  a  function 
of  TATB  content  for  the  slow  heating  rate.  The  simulated  results  are  now  closer  to  the 
experimental  results,  although  the  effect  of  including  the  temperature  dependence  of 
the  specific  heat  and  thermal  conductivity  has  been  to  reduce  the  rate  of  temperature 
increase  with  increasing  amounts  of  TATB.  The  simple  model  (Figure  3)  showed  an 
increase  in  reasonable  agreement  with  experiment,  whereas  Figure  15  shows  that  the 
improved  model  predicts  that  the  surface  temperature  is  almost  independent  of  TATB 
content.  Figure  16  shows  time  to  reaction  as  a  function  of  TATB  content  at  the  slow 
heating  rate.  The  simple  model  (Figure  4)  predicted  an  increase  in  time  to  reaction  as 
the  TATB  content  increased,  while  experimentally  there  is  a  decrease.  In  Figure  16 
the  improved  model  predicts  that  the  time  to  reaction  will  be  independent  of  TATB 
content,  and  agreement  with  experiment  is  considerably  better.  The  very  low 
experimental  value  at  20%  TATB  level  is  also  rather  suspect,  and  if  this  value  were 
higher  then  there  would  be  excellent  agreement  with  experiment. 

At  the  fast  heating  rate  the  improved  model  has  had  little  effect  on  the  simulated 
results,  although  the  small  changes  which  have  been  made  appear  to  have  shifted  the 
results  slightly  away  from  experiment.  Figure  17  shows  the  surface  temperature  at 
reaction  as  a  function  of  TATB  content.  The  simple  model  (Figure  7)  showed  a  slight 
increase  in  temperature  as  TATB  content  increased.  Figure  17  however  shows  that  the 
improved  model  predicts  a  very  slight  decrease  with  TATB  content.  A  similar 
discrepancy  is  shown  in  Figure  18,  which  shows  time  to  reaction  as  a  function  of  TATB 
content.  The  simple  model  (Figure  8)  showed  a  slight  increase  in  time  to  reaction  as 
the  TATB  content  increased,  while  experiment  shows  an  even  greater  increase.  Figure 
18  however  shows  that  the  improved  model  predicts  a  very  slight  decrease  in  time  to 
reaction  as  TATB  content  is  increased.  One  possible  explanation  for  this  discrepancy 
was  thought  to  be  the  use  of  equation  (9),  which  may  not  be  valid  at  the  higher  levels 
of  TATB  concentration.  Hence  the  simulated  results  shown  in  Figures  17  and  18  were 
recalculated  assuming  that  the  thermal  conductivity  of  the  mixture  was  given  by  the 
mass  average  of  the  values  for  the  individual  explosives,  as  used  by  McGuire  and 
Tarver  [3].  This  altered  the  simulated  results  by  less  than  0.5%  however.  We  currently 
have  no  explanation  for  the  failure  of  the  improved  model  to  predict  better  agreement 
with  experiment  for  TATB  at  the  fast  heating  rate. 
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Figure  11:  Surface  temperature  at  reaction  as  a  function  of  PETN  content.  The  open  circles 
are  the  experimental  points  and  the  filled  circles  are  the  predictions  of  the  improved  model. 


Figure  12:  Time  to  ignition  as  a  function  of  PETN  content.  The  open  circles  are  the 
experimental  points  and  the  filled  circles  are  the  predictions  of  the  improved  model. 
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SSCB  Tests  of  RDX/PETN/Elvax  210 


Figure  13:  Surface  temperature  at  reaction  as  a  function  of  PETN  content.  The  open  circles 
are  experimental  points  and  the  filled  circles  are  the  predictions  of  the  improved  model. 


Figure  14:  Tune  to  reaction  as  a  function  of  PETN  content.  The  open  circles  are  the 
experimental  points  and  the  filled  circles  are  the  predictions  of  the  improved  model. 
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Figure  15:  Surface  temperature  at  reaction  as  a  function  ofTATB  content.  The  open  circles 
are  the  experimental  points  and  the  filled  circles  are  the  predictions  of  the  simple  model. 
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Figure  16:  Time  to  reaction  as  a  function  of  TATB  content.  The  open  circles  are  the 
experimental  points  and  the  filled  circles  are  the  predictions  of  the  improved  model. 
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Figure  17:  Surface  temperature  at  reaction  as  a  function  ofTATB  content.  The  open  circles 
are  the  experimental  points  and  the  filled  circles  are  the  predictions  of  the  improved  model. 


Figure  18:  Time  to  reaction  as  a  function  ofTATB  content.  The  open  circles  are  experimental 
points  and  the  filled  circles  are  the  predictions  of  the  improved  model. 
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6.  Conclusion 


We  have  described  the  development  of  a  one-dimensional  model  designed  to  simulate 
heat  flow  in  a  two-component  explosive  mixture.  The  model  was  implemented  in  two 
stages,  and  at  each  stage  the  results  of  the  simulations  were  compared  with 
experimental  values  obtained  from  SSCB  tests.  In  the  first  stage  we  neglected  the 
temperature  dependence  of  the  thermophysical  constants  of  the  explosives,  used  a 
zero  order  kinetic  scheme  for  their  thermal  decomposition,  and  only  considered 
melting  of  the  RDX  component.  Even  with  these  limitations,  the  simulated  results 
were  found  to  reproduce  the  overall  trends  shown  by  experiment,  the  only  exception 
being  the  time  to  reaction  for  RDX/TATB/Elvax  210  at  the  slow  heating  rate,  where 
the  experimental  results  decreased  slightly  as  the  TATB  content  increased,  while  the 
simulated  results  predicted  an  increase  in  time  to  reaction  as  TATB  content  increased. 

In  the  second  stage  we  added  various  refinements  to  the  model.  We  included  the 
temperature  dependence  of  the  thermal  conductivity  and  specific  heat  of  RDX  and 
TATB,  and  a  term  to  model  melting  of  the  PETN.  We  also  used  a  first  order  kinetic 
scheme  for  the  decomposition  of  each  of  the  explosives.  In  general,  we  found  that 
these  refinements  greatly  improved  agreement  with  the  experimental  results.  For 
RDX /TATB /Elavax  210  at  the  fast  heating  rate  however  the  improved  model  resulted 
in  a  very  slight  change  in  the  simulated  results  away  from  experiment. 

Our  improved  model  also  considered  the  effect  of  a  multi-term  kinetic  scheme  to 
describe  the  decomposition  of  the  RDX,  and  this  was  found  to  provide  poor 
agreement  with  the  current  SSCB  results.  As  discussed  in  section  4,  this  is  believed  to 
be  related  to  the  differing  amounts  of  confinement  in  the  SSCB  and  ODTX  tests. 

In  conclusion,  our  calculations  clearly  illustrate  the  need  to  include  the  temperature 
dependence  of  the  thermal  properties  of  the  material,  and  a  kinetic  decomposition 
scheme  appropriate  to  the  degree  of  confinement,  before  good  agreement  between 
simulated  and  experimental  results  can  be  obtained. 
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